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We study the exclusive photoproduction of a lepton pair off the proton with the aim of studying the 
proton quark structure via the Generalized Parton Distributions (GPD) formalism. After deriving 
the amplitudes of the processes contributing to the 7 P —> P'e + e ~, the Timelike Compton Scattering 
and the Bethe-Heitler process, we calculate all unpolarized, single- and double- beam-target spin 
observables in the valence region in terms of GPDs. 

PACS numbers: 


I. INTRODUCTION 

The scattering of light on matter, which can generically 
be called Compton scattering, is a powerful tool to inves¬ 
tigate its inner structure. Nowadays, understanding the 
structure of hadrons in terms of quark and gluon (par- 
tons) degrees of freedom, i.e. the basic constituents of 
matter known to this day, is the subject of an intense re¬ 
search effort. Only these past fifteen years or so, thanks 
to the emergence of high intensity, high energy (multi- 
GeV) and high duty-cycle lepton accelerators, Compton 
scattering at the partonic level starts to be investigated 
experimentally in an efficient way. 

A particular case of Compton scattering at the par- 
tonic level is the Deeply Virtual Compton Scattering 
(DVCS) process on the proton P, i.e. y*P -T 7 P' where 
the initial virtual photon 7 * is radiated from an incom¬ 
ing lepton beam (see Fig. [l]-top). It is of particular in¬ 
terest as the amplitude of the process allows to access 
some essential operators of Quantum Chromo-Dynamics 
(QCD, the fundamental theory governing the interactions 
of quarks and gluons). Indeed, at sufficiently large vir¬ 
tuality of the initial photon (Q 2 = (k — k') 2 ), the DVCS 
amplitude can be factorized into an elementary “hard” 
(perturbative) scattering process 7 *g -A 7 # (where q is a 
quark of the proton), exactly calculable from Quantum 
Electro-Dynamics as well as perturbative QCD, and a 
“soft” (non-perturbative) QCD bilocal matrix elements. 
The Fourier transforms into momentum space of these 
QCD matrix elements are the so-called Generalized Par- 
ton Distributions (GPDs). In DVCS on the nucleon, 
at QCD leading-twist order there are four quark helic- 
ity conserving GPDs (H , E, H and E) which can be 
accessed and which correspond to the four independent 
helicity-spin transitions between the initial quark-proton 
system and the final one. The GPDs contain a wealth 
of information on the partonic structure of the proton: 
the longitudinal momentum and transverse space distri¬ 
butions of the quarks and gluons, the correlation between 
these momentum and space distributions, sensitivity to 
the quark-antiquark content in the proton, the quark or¬ 
bital momentum contribution to the proton spin, etc... 
We refer the reader to the reviews m on GPDs for the 


details of the formalism and of their properties. 




Figure 1: Top: QCD leading-twist DVCS diagram; Bottom: 
QCD leading-twist TCS diagram (the crossed diagrams of 
both processes are not shown). The dashed line illustrates 
the factorization between 


DVCS is currently widely investigated, experimentally 
as well as theoretically. Ref. [4] compiles all existing data 
in the valence region and presents some first informa¬ 
tions on the partonic structure of the proton that can 
be extracted, within some approximations, from the first 






2 


DVCS data through the GPD formalism. One example 
of a pioneering result is the first quantitative evaluation 
of the increase of the transverse size of the proton as 
smaller longitudinal momentum fractions of partons are 
probed. 

However, the GPD information is difficult to extract: 
DVCS cross sections are small (of the order of pb), there 
is the competing Bethe-Heitler (BH) process which leads 
to the same final state ePy but where the final photon is 
emitted from the incoming or scattered lepton and which 
is therefore not related to GPDs but nevertheless inter¬ 
feres. Furthermore, there is the need of measuring a se¬ 
ries of spin (beam and/or target) observables to constrain 
the different GPDs, etc... It would therefore be useful to 
investigate if supplementary and/or complementary con¬ 
straints on GPDs could be obtained from processes other 
than DVCS. In this spirit, we investigate in this arti¬ 
cle the related process of Timelike Compton Scattering 
(TCS) which corresponds to the exclusive photoproduc¬ 
tion of a lepton pair on the proton: 7 P —>• P'y* ^ e + e - 
and which is displayed in Fig. [l]-bottom. Like for DVCS, 
at sufficiently large virtuality of the final virtual photon 
(Q' 2 = (k + k') 2 ), it is predicted that the process factor¬ 
izes and is sensitive to GPDs, the same ones accessed in 
DVCS. 

The TCS process was originally investigated in terms 
of GPDs about ten years ago in Ref. [5]. In this pio¬ 
neering work, analytical formulas in terms of GPDs were 
derived for the unpolarized and the circularly polarized 
beam cross sections of the process 7 P —)> P'e + e _ , i.e. on 
a proton target. Very recently, a second article continued 
the investigation by studying the linearly polarized beam 
cross section 0 . 

However, in order to obtain simple analytical expres¬ 
sions, a few approximations were used in the calculation 
of the TCS amplitude (for instance mass correction terms 
of the order of m 2 /Q' 2 where m is the mass of the proton 
were neglected). In the present work, we waive some of 
these approximations and present calculations of differ¬ 
ent observables. In addition, besides unpolarized cross 
sections, we study all single and double beam-target spin 
observables. We focus in this article on a proton target. 

This article is organized as follows: in the next sec¬ 
tion, we present the general theoretical formalism of the 
TCS process, in particular the expression of the QCD 
leading-twist amplitude in terms of GPDs, and of the ac¬ 
companying BH process. We discuss some experimental 
considerations in the third section. In the fourth sec¬ 
tion, we present our numerical results for the unpolarized 
cross section of the 7 P —* P'e + e~ process and we com¬ 
pare them to the previous work of Ref. [51. In the fifth 
section, we present our results for single spin observables 
(beam and target) and we compare our work for the beam 
polarization observables to the results of Ref. [5]. In the 
sixth section, we present our results for the double-spin 
beam-target observables. For each case, we will show the 
dependence of the observables on different GPDs and its 
sensitivity on different kinematics. In the seventh section 


we show the impact of next-to-leading-twist corrections 
on the cross sections and on the asymmetries. We will 
conclude in the eighth section. 


II. FORMALISM 

We are studying the process: 

7 («)P(p)-> P\p')e-(k)e + (k') (1) 

in a GPD framework, i.e. when the final photons virtu¬ 
ality Q' 2 = (k + k') 2 is sufficiently large and the proton 
momentum transfer t = (p — p') 2 is sufficiently small so 
that the factorization illustrated in Fig. [l]-bottom can be 
applied. From DVCS, Deep Inelastic Scattering (DIS) 
and Drell-Yan analysis and experiments, it is believed 
that Q' 2 > 2 GeV 2 and —t< 1 GeV 2 (or < 30%) 
should define such a reasonable phase space. Regarding 
Q /2 , one should also avoid regions where one can have the 
production of vector mesons, decaying into e + e - pairs 
(for instance, the broad p'(1700)). Finally, one should 
consider the squared center-of-mass energy of the incom¬ 
ing photon and target proton s = (q + p) 2 > 4 GeV 2 , in 
order to minimize possible contributions from the Dalitz 
decay of proton resonances. 


A. Kinematics 

We will use the notation of Ji [7] for GPDs, i.e. GPDs 
depend on the three variables x, £ and t where the quark 
longitudinal momentum fractions x and £ are defined 
w.r.t. the average proton momentum P and proton mo¬ 
mentum transfer A, respectively. We therefore define: 


p = \(p + p'), 

(2) 

A = {p' ~p) = (q~ q'), 

(3) 


and we also introduce the average photon momentum 

Q=\{q + j)- ( 4 ) 

GPDs are matrix elements of QCD operators which are 
defined at equal light-cone time. It is therefore conve¬ 
nient to use a frame where the q and P momenta are 
collinear along the z- axis and in opposite directions. We 
define the lightlike vectors along the positive and nega¬ 
tive z directions as: 

r = P + /V2(l, 0 , 0 , 1 ), (5) 

= 1/P + ■ l/v / 2(l, 0,0, —1), ( 6 ) 

and we define the light-cone components a ± by a ± = 
(a 0 dz a 3 )/ y/2. We have p 2 = n 2 = 0 and p • n = 1. In this 
light-cone frame, introducing the + components of the 
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A and g four-vectors, £ and £' respectively, the various 
four-vectors involved can be decomposed as: 

_ o 
777 

P»=pn + n M , (7) 

= -i'r -^, n \ ( 8 ) 

A M = -2£p»* + |mV + A^, (9) 

with to 2 = m 2 — and to is the proton mass. We 
relate the final photon virtuality to the average photon 
momentum: 



q 2 = Q' 2 /2 - A 2 /4. (10) 


We have to relate the light-cone momentum fractions 
£ and £' to the kinematical variables which are experi¬ 
mentally accessible. To do so, we introduce the variables 
£ and £': 


r = - 


2P.q 




A .q 
2 P.q 


-Q' 2 + A 2 /2 
2(s — to 2 ) + A 2 — Q' 2 ’ 

Q ' 2 

2(s — to 2 ) + A 2 — Q' 2 ’ 


( 11 ) 

( 12 ) 



The light-cone momentum fractions £ and £' are re¬ 
lated to the kinematical variables £ and £' by: 


c = c ■ 




1 + £ /2 m 2 /q 2 
1 — £ /2 fh 2 /q 2 ’ 
2 


1 + \/l — 4£ /2 m 2 /q 2 


(13) 

(14) 


In the asymptotic limit, where mass and A terms are 
neglected relatively to Q /2 , we have: 


l=e = -f = -e' = 


Q 


,/2 


2 s - Q* 


(15) 


B. Timelike Compton amplitude 


The two diagrams of Fig. [2] have to be calculated. The 
leading order amplitude reads: 
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Figure 2: Leading order diagrams for the process yg —>> yg. 
We indicate between parenthesis the + components of the 
four-vectors. 


with, in the Bjorken limit where £ = — £', 


H 


TCS 




(18) 


= 1 (-g^u)± f dx ( -L_ + 1 ) 

2 ^ J \x — £ — ie x + £ + 2 ey 

-l 
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The full TCS amplitude, corresponding to the diagram 
of Fig. [l]-bottom (plus the associated crossed diagram) 
reads then 


(16) 


~ ~ /\ 
H{x,i,t)u{jp')^ u(p) + E(x,Z,t)u(p')i 5 — u(p) ) , 


where we used the metric: 

( 9/iv)I. = -9\iv + V»n v + p u n^ , (19) 

(^T^)_L ^u/ja/3 pP. 

The GPDs entering Eq. [18] are proton GPDs, i.e. they 
read, in terms of quark flavors: 

Htcs(x,U) = \a ulv + \n dlv + \n slv - ( 20 ) 


jiTCS 


—n u ( k ) T V T) H i 


TCS 

/JLV 5 


(17) 


In this work, we will take the GPD parametrizations 
from the VGG model B IMS, which are summa- 
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rized in Ref. [4 and based on the Radyushkin double¬ 
distribution ansatz for the (x,£)-dependence [TTUT3] and 
on a Reggeized ansatz for the ^-distribution [lj TO] , At 
a couple of instances, in order to estimate the model de¬ 
pendence of our calculations, we will use the factorized 
ansatz for the ^-dependence of the H GPD [ 8 ]. Also, we 
will occasionally study the sensitivity of observables to 
the so-called D-term m, which is included in the VGG 
model and whose parametrization can be found as well 
in Refs. ni iMo]. 


C. Gauge invariance 



The TCS amplitude is not exactly gauge invariant, 
only in the Bjorken limit. To restore gauge invariance, 
twist-3 corrections of order A±/Q are needed. As a first 
step in adressing this issue and estimating its effect, we 
propose to add a correction term to the twist -2 vector 
part tensor as follows: 


H^ = H^ 0 


pn 

2P~ 


= ' (Aj 


+ 


2 Pq 

pH pu 

'±{P-q) 2 


• (A ± ) a • H» x 


LO 


(AjJ* • (A_l)a 


TJKV 

■ h lo 


■ h kX 


LO 5 


( 21 ) 



where stands for the tensor of Eq. 18 This is a 

generalization of the prescription proposed in Refs. p[T5] 
for DVCS (we also refer the reader to Refs. fTSHTS] for 
further discussions on the issue of gauge invariance in 
the DVCS amplitude). One can readily check that 
respects gauge invariance both w.r.t. initial and final 
photons, i.e. q /1 H fJ/l/ = 0, and = 0. The impact of 

adding this correction to the observables is shown at the 


end of this paper, in the section VII 


Figure 3: The Bethe-Heitler diagrams. 


can be considered to be known with good accuracy. 
In this work, we take the parametrizations issued from 
Refs. 19> 20 ] . 


D. The Bethe-Heitler amplitude 


The TCS process is accompanied by the BH process, 
involving the two diagrams which are presented in Fig. [3] 
Their amplitude reads: 

T BH = -^u(p')r'u(p)e»(q) ( 22 ) 

m v{k%(n) 

with the virtual photon-proton electromagnetic vertex 

matrix 

in v P A 

r l ' = 'fF 1 {t) + ——e-F 2 (t). (24) 

2 m 

The BH amplitude depends on the proton Dirac and 
Pauli form factors F\ (t) and which, at small t 


E. Cross section 

At fixed beam energy E 1 or longitudinal momentum 
transfer £, there are four independent kinematical vari¬ 
ables for the process 7 (q)P(p) —>• P'(p')e~ (k)e+ (k'). A 
natural choice is to take: Q' 2 and t = A 2 that we already 
defined, and the two angles 6 and 0 of the electron in the 
7 * center-of-mass (with the V-axis along the direction 
of the 7 * in the 7 * — P' center of mass). We illustrate 
in Fig. [4] the different variables involved. In addition, we 
display the polarization angle T between the polarization 
vector of the incoming photon and the scattering plane 
in the 7 * — p C.M.. 

The 4-differential unpolarized cross section is then ex- 
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Figure 4: Top panel: scheme of the reaction plane 7 P —>> 7 *P' 
(for TCS) in the 7 P CM frame and illustration of some of the 
kinematic variables. Red arrows show the polarization vectors 
of the photon beam and proton target. 

Bottom panel: decay angles in the e _ e + CM frame. 0 and 6 
are respectively the azimuthal and polar angles between the 
e~ direction and the virtual photon direction in the 7 * — P' 
CM frame. 


pressed as: 


d A c 


( 7 P -T P'e + e~) 


dQ ,2 dtdPt 

11 1 

(27t) 4 64 (2mP 7 ) 2 


(25) 


rBH 


r TCS 12 


where | T BH -\-T tcs | 2 is averaged over the target proton 
helicities and beam polarizations and summed over the 
final proton helicities. 


III. EXPERIMENTAL CONSIDERATIONS 

The only data for the process 7 P —>> P'e~e + which are 
in the phase space of concern for our study, have been 
collected and analyzed a few years ago by the CLAS col¬ 
laboration using the ~ 6 GeV electron beam of Jeffer¬ 
son Lab. The data were actually obtained in “quasi¬ 
photoproduction” mode. This means that the scattered 
electron from the beam is not detected in CLAS and 
is considered to be in almost the same direction as the 
beamline. This results in very low Q 2 electroproduc¬ 
tion, i.e. “quasi-photoproduction”. This pilot analysis 
can be found in Ref. [ 21 ]. The « 6 GeV beam energy, 
combined with the ~ 10 34 cm - 2 s -1 luminosity, allowed to 
reach maximum Q' 2 values of 3 GeV 2 , which corresponds 
to an invariant mass of the e + e - system M e + e - of ~ 1.8 
GeV. This is close to the mass of several vector mesons 
decaying into e + e - , in particular the broad p'(lTOO). In 
order to have a TCS interpretation as clean as possible, it 
is of course advisable to avoid such resonances which con¬ 
tribute to the process 7 P —>> P'e~e + . The data analysis 
of Ref. [ 21 ] is therefore difficult to interpret in terms of 
GPDs but it nevertheless demonstrates that the process 
7 P P'e~e + can be measured at JLab. 

The forthcoming JLab energy upgrade to 11 GeV al¬ 
lows to explore a Q 12 region between 4 and 9 GeV 2 which 
corresponds to a M e + e - region between 2 and 3 GeV, i.e. 
a vector meson resonance-free region between the p' and 
the J/\h. In the case of CLAS 12 , the upgraded CLAS 
detector associated to the JLab energy increase, there is 
also a luminosity gain of a factor 10. These upgrades 
have led to the first dedicated TCS accepted proposal at 
JLAB [ 22 ]. It will use a similar “quasi-photoproduction” 
technique as used in the pioneer 6 GeV analysis. In ad¬ 
dition, there will be the improvement of the detection 
of the low Q 2 scattered electron via a dedicated tagging 
equipement, supplementing CLAS 12 . This allows, be¬ 
sides the measurement of unpolarized cross sections, to 
obtain linearly polarized photons observables, by measur¬ 
ing the azimuthal dependence of the scattered electron. 
Finally, with a polarized electron beam, which is avail¬ 
able at JLab, one can access circularly polarized photons 
observables. It is therefore expected that in the next few 
years, numerous 7 P —)► P'e~e + data which can lend to 
GPD interpretation will be available. 

We show in Fig. [5] the kinematical domain which can 
be accessed with the upgraded JLab. We display in blue 
the (£, Q /2 ) phase space accessible for TCS with an 11 
GeV electron beam, assuming that the real photon is 
provided by bremsstrahlung of the electron and that its 
energy is in E 1 G [5, 11] GeV. We have applied two cuts: 
Q' 2 G [4,9] GeV 2 and —t G [0,1] GeV 2 . The motivations 
are respectively to stay in the region free of vector mesons 
resonances and minimize higher twist corrections to the 
TCS formalism, which grow with We overlap in red 
in this same figure the (£, Q 2 ) phase space accessible with 
an 11 GeV beam for DVCS. We have applied the cuts: 
—t G [0,1] GeV 2 , for the same reason as for TCS, and 
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s > 4 GeV 2 in order to stay above the baryon resonance 
region. 

One notes the large intersection between the DVCS 
and the TCS phase spaces. Measurements of observables 
sensitive to GPDs in the common (£, Q 2 ) region by both 
processes should bring strong constraints on the extrac¬ 
tion of GPDs and tests of factorization and universality. 



Figure 5: Kinematical domain accessible as a function of £ and 
Q 2 for DVCS (red plain surface) and in £ and Q' 2 for TCS 
(blue dotted surface) with an 11 GeV electron beam. For 
DVCS, the cuts —t G [0,1] GeV 2 and s > 4 GeV 2 have been 
applied and for TCS, the cuts E 1 G [5,11] GeV, — t G [0,1] 
GeV 2 and Q' 2 G [4,9] GeV 2 have been applied. 


We now present our results for the calculations of the 
unpolarized cross sections, single spin asymmetries and 
double-spin asymmetries, respectively in sections |IV[ |V| 
and |VI| In these sections, all calculations will be done 
in the Bjorken limit and at leading order (we refer the 
reader to Refs. [23H25] for works on Next-To-Leading or¬ 
der corrections in a s to the TCS amplitude). In section 
m we will study the effect of keeping the exact kinemat¬ 
ics presented in section |II A| and of the gauge-invariance 


restoration prescription described in section IIC 


IV. UNPOLARIZED CROSS SECTION 

We discuss in this section the unpolarized cross section 
of the yP —»• P'e + e~ process, which therefore includes 
the BH and TCS processes. Fig. [6] shows our calculation 
of the (^-dependence of the 4-fold differential cross sec- 
^° n dQ ' 2 dt d^d^cosQ ) £ = _ ^ = ^ GeV 2 , Q 2 = 7 

GeV 2 and for different 0 values. The 0-shape is strongly 
dependent on the 0 value. As 0 tends to 0°, the 0 dis¬ 
tribution peaks towards 0=180° and as 0 tends to 180°, 
the 0-distribution peaks towards 0=0° (or 360°). There 
is a smooth transition between these two singularities for 


the intermediate 0 values. For instance, at 0=90°, the 
calculation shows only two small “bumps” at 0=0° and 
0=180°. 

These particular shapes are due to the BH process and 
its singularities. Indeed, in the BH diagrams of Fig. [3| 
when the electron (positron) is emitted in the direction 
of the initial photon, i.e. 0=0° (0=180°), the propa¬ 
gator of the positron (electron) becomes singular and 
creates a peak in the 0 distribution at 0=180° (0=0°). 
Intuitively, 0=0° (0=180°) forces all particles to be in 
the same plane, i.e. 0=180° (0=0°). The kinematics 
0=0°, i.e. the electron is in the direction of the photon 
beam, corresponds to 0=180° because the virtual photon 
is emitted by the positron, not the electron (see Fig. [6|. 

We display also in Fig.[6]the contribution of TCS alone. 
In this calculation, we have used only the GPD H. The 
inclusion of the other GPDs barely changes the curves. 
In contrast to the BH, the TCS is almost flat in 0 for 
all 0 values. It is clear that the process yP -T P'e + e~ 
is largely dominated by the BH. There is never less than 
an order of magnitude between BH and TCS. 

In Fig. [6j we show the curves BH+TCS for a series of 
0 angles as well as the BH alone for 0=45° and 0=90°. 
Only at 0=90°, where one is far from the two BH sin¬ 
gularities, we have a visible difference between the two 
curves and therefore a sensitivity to TCS. It is of the or¬ 
der of 30% at 0=180°. As one gets closer to one of the two 
BH singularities (0=45° for instance), the two curves BH 
and BH+TCS are essentially indistinguishable and there 
is no sensitivity to TCS. 

Finally, we show in Fig. [6j our calculations of BH+TCS 
(and of BH alone) for 0 integrated over the range 
[7r/4, 37 t/ 4]. In order to maximize count rates, from an 
experimental point of view, it is interesting to integrate 
over 0. We still have a sensitivity to TCS, however it is 
of the order of 5%, i.e. lesser than at fixed 0=90°: the 
integration over 0 dilutes the sensitivity to TCS. 

We display in the top panel of Fig. [ 7 ] the (In¬ 
dependence, at £ = 0.2 and —t = 0.4 GeV 2 , of d<j/dQ r2 dt 
for the BH and TCS processes. The 2-fold cross section 
has been integrated over the decay angles: 0 G [0, 27r] 
and 0 G [J, ^]. We also display in this figure the results 
of the analytical formulaes of Ref. [5]. For TCS, we have 
of course used the same GPD parameterization for both 
calculations (only H in this case). 

In order to better appreciate the comparison, we dis¬ 
play in the bottom panel of Fig. [7| the ratios of our calcu¬ 
lations to those of Ref. [5]. The agreement is excellent for 
the BH process, as it should, since there was no approx¬ 
imation done in the derivation of the analytical formula 
for this process in Ref. [5]. For TCS, we show two curves 
for the ratio. Indeed, in Ref. [5], in the analytical formula 
for the TCS process, the nucleon mass was neglected in 
the phase space factor, i.e. the cross section is propor¬ 
tional to -4 rather than — ^2 (this is not the case for 
the BH process where the phase space factor is exact 
in Ref. [5]). In the bottom panel of Fig. [7| we plot the 
TCS cross section of Ref. |5j with (blue dotted curve) and 











7 


- bh+tcs, 0 e(jr/4,3ji/4) 

bh+tcs, 0=10° 
bh+tcs, 0=45° 
bh+tcs, 0=90° 
bh+tcs ,0=135° 

- bh+tcs, 0=170° 

bh, 0 <E(ji/4,3ti/4) 
bh, 0=45° 
bh, 0=90° 

- tcs, 0 <E(ji/4,3jt/4) 

tcs, 0=10° 
tcs, 0=45° 
tcs, 0=90° 

- tcs, 0=135° 

tcs, 0=170° 



process 7 P P'e + e~ (for BH + TCS and TCS alone) as 

a function of £>, at £= 0 . 2 , —t — 0.4 GeV 2 , Q' 2 = 7 GeV 2 
for different fixed 6 values: 10°, 45°, 90°, 135°, 170° and for 
6 integrated over [7r/4, 37 t/4]. For two values of 0 (45° and 
90°) and for 0 integrated over [7r/4, 37 t/4], we also show the 
calculation of the BH alone. 


without (red dashed curve) the nucleon mass term in the 
phase space factor. In this way, one can distinguish the 
differences between the two calculations originating from 
the trivial phase space factor from those, more subtle, 
coming from the TCS amplitude. It is seen in this figure 
that at the lowest Q' 2 values and at the presently consid¬ 
ered kinematics, the difference in the cross section calcu¬ 
lations depending on the prescription taken for the phase 
space factor can reach more than 10%. In both cases, it 
is also seen that the difference between the Berger et al.’s 
calculations and ours diminishes as Q' 2 increases, as ex- 

2 

pected since terms of the order of ^2 were dropped in 


the TCS analytical formula of Ref. [5]. 



J_ 1 1 1 _ 1 I 1 1 1 1 I 1 1 _ 1 _,_l_ 1 1 1 1 I 1 1 1 1 _L 

4 5 6 7 8 9 

Q ’ 2 (GeV 2 ) 



O ’ 2 (GeV 2 ) 


Figure 7: Top panel: the two-fold differential cross section 
da/dQ' 2 dt for BH (highest curves) and TCS (lowest curves) as 
a function of Q' 2 for £ = 0.2 and —t = 0.4 GeV 2 . The calcu¬ 
lations have been integrated over £> G [0, 2tt] and 0 G [f, ^]. 
The solid curves are our calculations while the dashed curves 
are from Ref. 5l. Bottom panel: ratio of our calculations to 
those of Ref. o . For the TCS process, the blue dotted curve 
contains the nucleon mass term in the phase space factor in 
the calculation of Ref. [5] while the red dashed curve does not. 


To end this section concerning the unpolarized cross 
section, we show in Fig. [8] the influence of the D-term 
on the three-fold differential cross section da/dQ' 2 dtdct>, 
calculated for £ = 0.2, —t = 0.4 GeV 2 and with 6 in¬ 
tegrated over [f, It modifies the amplitude of the 
cross section at 0 = 0° and 0 = 180° by about 10%. 



































<|> (deg.) 


Figure 8: The three-fold differential cross section 

da / dQ' 2 dtdfi as a function of 4> with (red dashed curve) and 
without (black solid curve) the D-term, calculated at £ = 0.2, 
—t = 0.4 GeV 2 and with 0 integrated over [^, ^f-\. 


V. SINGLE SPIN ASYMMETRIES 

We now turn to the single spin asymmetries: beam 
or target. Photons beams can be polarized linearly or 
circularly. Target polarization vectors can be oriented 
along the x— , y— or z—axis in the 7 P -+ y*P' plane 
(see Fig. [4|. For polarization observables, we will calcu¬ 
late spin asymmetries and, following notations used for 
DVCS, we will tag them i.e. with two indices i and 
j. The first index i refers to the polarization type of the 
beam: U for an unpolarized beam, © for a circularly po¬ 
larized beam and £ for a linearly polarized beam. The 
second index j refers to the polarization of the target 
and can take the values J7, x, y or z, with obvious mean¬ 
ings. In this section, dedicated to single-polarization ob¬ 
servables, we will therefore consider successively the five 
independent asymmetries Agjj, Aqu, Ajj x , Ajj y and Ajj z - 


A. Linearly polarized photons 


We introduce the angle 4/ between the polarization vec¬ 
tor of the photon and the plane spanned by the photon 
beam and the (e + e - ) system, which contains the x-axis 
(see Fig. [4|. Then, we define: 


Aiu{^) 


W - <Ty(*) 

o x (*)+° y (*r 


(26) 


where cr x (a y ) stands for the 4-fold differential cross sec- 

tions dQ '2 dt d$d(cos8)dV with a P hoton beam P olarized 

in the x-(?/-)direction. We display in Fig. [9] the T- 


dependence of Am for £= 0 . 2 , —t = 0.2 GeV 2 , Q' 2 = 7 
GeV 2 for 0 = 45°, 0 = 90° and 0 integrated over 
[7r/4,37r/4] and 0 = 0° and = 10°. The approximate 
shape of the asymmetry is a cos( 2 \h) which is reminiscent 
of the modulation which is predicted for the so-called 
E asymmetry in single meson photoproduction. This 
cos( 2 T) modulation appears explicitely in the analytical 
expressions of Ref. [ 6 ,. 

We note that the BH alone produces an asymmetry. 
It is actually the dominant contribution. The TCS pro¬ 
duces only variations of the amplitude at the percent level 
around the BH, making this observable not very favor¬ 
able to study TCS and GPDs. The amplitude and phase 
of the asymmetry depend strongly on the decay angles: 
it is the strongest as 0 approaches 90° and the phase in¬ 
creases as (p increases. This phase shift due to </> is also 
apparent in Ref. j6j. 

In Fig. [9j we have used only the contribution of the H 
GPD for TCS. In Fig. [lOj we show the ©dependence of 
the Aw asymmetry for BH alone, BH+TCS (with only 
H ) and BH+TCS (with H+H). Calculations have been 
done for £= 0 . 2 , Q' 2 = 7 GeV 2 , 0 = 0 °, 4/ = 180° and 
for 6 = 45°, 6 = 90° and 0 integrated over [7t/4, 37t/4]. 
Depending on t, we notice some sensitivity of the A^u 
observable to the GPDs H and H. 

In Fig. [TTJ we show the peculiar ^-dependence of Am 
at 4/ = 0° at our standard kinematics. We also show the 
(small) influence of the D-term. 


- BH, 0=90°, 4)=0° — BH, 45°<0<135°, (j)=0 o 

BH+TCS, 0=90°, (p=0° BH+TCS, 45°<()<135°, <p=0° 


BH, 0=45°, (j)=0° . BH, 0=90°, 4>=10° 

BH+TCS, 0=45°, (J)=0 o 



T (deg.) 


Figure 9: Am asymmetry as a function of T for BH and for 
BH+TCS (only GPD H contribution) at £=0.2, —t — 0.4 
GeV 2 , Q /2 = 7 GeV 2 and for different sets of 6 and £> angles. 
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- 6=90°: BH 

BH+TCS, H 
BH+TCS, H+H 
0=45°: BH 
BH+TCS, H 
BH+TCS, H+H 
45°<0<135 O : BH 
BH+TCS, H 
BH+TCS, H+H 



Figure 10: Am asymmetry as a function of t for £=0.2, Q /2 « 
7 GeV 2 , T = o°, £> = 0° and <9 = 45°, 6 = 90° and 0 e 
[7t/4, 37t/4]. 


B. Circularly polarized photons 


We define: 


Aqu = 


— & 


cr+ + a 


(27) 


where a± stand for the 4-fold differential cross sections 
d Q f2 dt ^ d(cos 6 ) f° r the two Photon spin states, right and 
left polarized. 

We display in Fig. fl2] (top panel) our results for Aqjj 
as a function of 0 at Q' 2 = 7 GeV 2 , £ = 0.2, —t = 0.4 
GeV 2 for 0=45°, 90° and 6 integrated over [45°, 135°]. 
The TCS is calculated here with only the H GPD. In all 
kinematics, we obtain a sin(0) shape with a significant 
amplitude, up to ~ 25%. We observe that the BH doesn’t 
produce any asymmetry. Any signal therefore reflects a 
contribution from TCS. This is due to the fact that, as 
was shown in Ref. |5], this observable is sensitive to the 
imaginary part of the amplitude and that the BH am¬ 
plitude is purely real. The amplitude of the asymmetry 
depends on 6. It is maximal for 6= 90° where BH is min¬ 
imal and it decreases as 0 tends to 0= 0° (or 180°). Since 
the BH does not produce on its own an asymmetry, one 



(deg.) 


Figure 11: Am asymmetry as a function of <p for £ = 0.2, 
Q' 2 = 7 GeV 2 , -t = 0.4 GeV 2 and 0 e [tt/4, 3tt/ 4]. The black 
solid curve corresponds to BH only, the red dashed curve to 
BH+TCS with only the GPD H contribution and the blue 
dotted curve to BH+TCS with the GPD H and the D-term 
contributions. 


sees that the integration over 0 does not strongly reduce 
the signal. Such integration allows to maximize count 
rates. 

In Fig. 12 (top panel), we also compare our results 
to those of Ref. [5]. In all cases, our calculations pro¬ 
duce amplitudes a few percents larger. This might be 
attributed to some mass correction terms to the TCS 
amplitude which are present in our calculation and not 
in Ref. [5] . 

Fig. [ 12 ] (bottom panel) shows the Aqjj asymmetry 
for 6 integrated over [45°, 135°] using different GPDs 
parametrizations for TCS. 

In Fig. 13 we show for £ = 0.2, Q' 2 = 7 GeV 2 , 0 = 90° 
and 0 integrated over [45°, 135°], the ^-dependence of 
Aqjj and its sensitivity to different GPDs. We notice 
that the magnitude of Aqu increases with \t\ and that 
there is a sensitivity of this observable to all four GPDs, 
especially at large | t |. We also display in this fig¬ 
ure our calculation with the factorized ansatz for the 
^-dependence of the H GPD in order to illustrate the 
model-dependence of our results. 


C. Polarized targets 


We define: 


Mji = 


<r+ + a 


(28) 


where a ± stands for the four-fold differential cross sec- 
tions d Q , 2 dt d ^ d ( cos Q^ for the two target spin orientations 
+ and — along the axis i = x, y or z. 
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<)> (deg.) 



<t» (deg.) 



|t| ((GeV 2 ) 


Figure 13: The Aqu asymmetry as a function of t for 

BH+TCS at £ = 0.2, Q' 2 = 7 GeV 2 , 0 = 90° and 0 inte¬ 
grated over [45°, 135°]. TCS is calculated with different GPD 
contributions. 


In this figure, TCS is calculated with different GPDs. 
Depending on the value of £, the three asymmetries are 
sensitive to the GPDs H, H and E in various proportions. 
We also display in this figure our calculation with the 
factorized ansatz for the t-dependence of the H GPD in 
order to illustrate the model-dependence of our results. 


Figure 12: Top panel: The Aqu asymmetry as a function 
of £> for 0 = 45°, 0 = 90° and 6 G [45°, 135°]. The TCS 
is calculated here with only the H GPD. We compare our 
calculations to those of Ref. [5]. Bottom panel: The Aqu 
for 9 G [45°, 135°] using differents GPDs parametrizations for 
TCS. The calculations are done for Q /2 = 7 GeV 2 , £ = 0.2, 
—t m 0.4 GeV 2 . 


We show in Fig.[M]our results for the ^-dependence of 
A Ux , A Uy and A Uz for Q' 2 = 7 GeV 2 , £ = 0.2, —t = 0.4 
GeV 2 for 6 integrated over [J, ^\. Like for Aqu, it is 
advantageous to integrate over 0 , in order to maximize 
count rates, since the signal is barely reduced. The TCS 
is calculated with different GPD contributions. We ob¬ 
serve sin^ or cos (j) shapes with amplitudes between 10 
and 15%. Like for Aqu , the BH doesn’t produce any 
asymmetry and any non-zero asymmetry directly reflects 
the strength of GPDs. 


We show in Fig. 15 the t-dependence of Au x , Au y and 
Auz at 0=90°, 0° and 90° respectively, for the kinematics 
£ = 0.2, Q' 2 = 7 GeV 2 and 0 integrated over [J, ^]. 
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Figure 14: The Ajj x (left panel), Ajj y (central panel) and Ajj z (right panel) asymmetries as a function of 0 for £ = 0.2, Q /2 = 7 
GeV 2 , —t = 0.4 GeV 2 and for 0 integrated over [^, ^f]. The TCS is calculated with different GPDs contributions. 
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Figure 15: The Ajj x (left panel), Au y (central panel) and Ajj z (right panel) asymmetries as a function of t, at 0=90°, 0° and 
90° respectively, and for £ = 0 .2, Q' 2 = 7 GeV 2 , —t = 0.4 GeV 2 and 6 integrated over [J, ^]. TCS is calculated with different 
GPD contributions. 


VI. DOUBLE SPIN ASYMMETRIES 


A. Circularly polarized photons and polarized 
target 


We define the double-spin asymmetries: 


_ (a ++ + a—) - (a + ~ + a~ + ) 
lj cr++ + a — + cr+ - + cr _ + 


(29) 


where cr ±=t stand for the four-fold differential cross sec¬ 
tions d Q , 2 dt ^ d( cos Q) f° r the two photon beam spin states 
+ and — (first index) and the two target spin orientations 
+ and — (second index) along the target polarization 
axis. The first index of A refers to the polarization na¬ 
ture of the photon beam (i = i for a linear polarization 
and 7 = 0 for a circular polarization) and the second in¬ 
dex refers to the axis polarization of the target j = x,y 
or z. We present and discuss in the two following sub¬ 
sections our results for Aqj and Aij. 


Fig. 16 shows our results for Aq x , A Qy and A@ z , from 
left to right, as a function of (f) at the kinematics £ = 0.2, 
—t = 0.4 GeV 2 , Q' 2 = 7 GeV 2 . The top row shows the 
result of our calculations for 9 integrated over [J, 
with different GPD contributions to the TCS amplitude. 
The bottom row shows the same observables with only 
the GPD H contribution for different 6 angle sets. 


One notes that the BH process alone produces asym¬ 
metries in all cases. The 0-shapes of the asymmetries are 
complex and very dependent on 6. Also, in contrast to 
the single spin asymmetries that we studied in the previ¬ 
ous section, the 0-shapes are also very dependent on the 
specific GPDs entering the TCS process. One notes in 
particular important sensitivities to the IT, iT, E GPDs 
as well as to the D-term. Furthermore, one also notes the 
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sensitivity of these double-polarization observables to the 
ansatz used for the ^-dependence. 


Fig. 17 shows the ^-dependence of the double spin 
asymmetries A© x , Aq v and A@ z at 0 = 0°, 0 = 90° and 
0 = 0° respectively, at the kinematics £ = 0.2, Q' 2 = 7 
GeV 2 and 6 integrated over [7r/4, 37 t/ 4] . We also display 
in this figure our calculation with the factorized ansatz 
for the t-dependence of the H GPD in order to illus¬ 
trate the model-dependence of our results. The change 
of sign for A Qy for the factorized ansatz compared to 
the Reggeized ansatz is in particular remarkable. As can 
be seen in the top panel of Fig. [l6j this is due to the 
fact that the factorized model crosses the “zero” line be¬ 
fore (j) = 90° while the Reggeized model crosses it after, 
thus producing respectively positive and negative A 02/ ’s 
VII. CORRECTIONS TO THE 
LEADING-TWIST CALCULATIONS 


We evaluate in this final section two types of higher- 
twist corr ectio ns: keeping the exact kinematics presented 
i-e. | £ | 7^1 £ It^I £' Ml £' I, and adding the 


in section 


IIA 


gauge correction term of section [TTC] to the TCS tensor. 


Fig. 19 shows the Q ' 2 -dependence and the In¬ 


dependence of the ratio of the 2-fold cross differential 
cross sections da/dtdQ' 2 for Bethe-Heitler and for TCS 
(integrated over the decay angles: (j) G [0, 27r] and 6 G 
[f, T 1 ]) ca l cu l a t e d with the gauge invariance restoration 
term (dot-dashed curve for TCS) and with the exact kine¬ 
matics (dotted curve for TCS) to the asymptotic limit 
(i.e. Bjorken limit) calculation that we have presented 
so far. The calculation has been carried out for £ = 0.1 
and at —t = 0.4 GeV 2 (left panel) or at Q' 2 = 7 GeV 2 
(right panel). One sees that these ratios tend to 1 as Q' 2 
increases and as \t\ decreases, as expected. The effects 
of these corrections in the domains of current interest, 
covered in this figure, are of the order of a few percents 
for TCS. The exact kinematic corrections are more im¬ 
portant (always less than 6% though) than the gauge 
invariance restoring corrections. 

We show in Fig.[l9]that the impact of these corrections 
on the asymmetries that we discussed earlier is of the 
order of 0.1% to 1% and don’t affect the conclusions 
that we drew earlier. 


VIII. CONCLUSION 

In this work, we have studied the 7 P —)> P'e + e~ pro¬ 
cess in terms of GPDs in a regime where one can expect 
to access them. We have presented our derivations of 
the TCS and BH amplitudes both contributing to the 
yP -> P'e + e _ process and calculated all unpolarized, 
single beam, single target and double beam-target spin 
observables. 

We show that, since the TCS process is lower by several 
factors in the unpolarized cross section compared to the 
BH process, it is judicious to measure spin asymmetries 


at (j) = 90°. 


B. Linearly polarized photons and polarized target 

Fig. [TB] shows our results for the double-spin asymme¬ 
tries At x , Ai y and Ai z , from left to right, as a function 
of 0 at the kinematics £ = 0.2, —t = 0.4 GeV 2 , Q' 2 = 7 
GeV 2 with 0 integrated over [J, ^]. 

One notes that the BH process alone produces a null 
asymmetry in all cases, making this double-spin asym¬ 
metry particularly favorable to study TCS and GPDs. 
Ai y is mostly sensitive to the GPD H contribution while 
Ai x and Ag z show a sensitivity to the GPD H as well. 

which reveal a more direct sensitivity to GPDs. In par¬ 
ticular, the BH process alone doesn’t produce any signal 
for the target single spin asymmetries, for the circularly 
polarized beam single spin asymmetries and for the lin¬ 
early polarized photons and polarized target double spin 
asymmetries. These observables are therefore particu¬ 
larly favorable to directly measure GPD strength. We 
have shown in general that the various single- and double¬ 
polarization observables that we have calculated show 
different sensitivities to the four GPDs, which should ul¬ 
timately allow to disentangle them with some adequate 
GPD fitting algorithms. We also provided first estima¬ 
tions of higher-twist contributions such as keeping the 
exact kinematics and including a gauge invariance restor¬ 
ing term. The effects are at the level of a few percent on 
cross sections and spin asymmetries. 

A rich new experimental TCS program can be envi¬ 
sioned with the forthcoming JLab energy upgrade, which 
would complement the DVCS program already approved 
in order to access GPDs. All of the TCS observables 
that we calculated in this work should be measurable 
and can serve as a basis for developping experimental 
proposals. This work might also find some applications 
at higher energies, like at hadron colliders, such as LHC 
and RHIC, in ultraperipheral collisions |26j or at the pro¬ 
jected electron-ion collider EIC m- 
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Figure 16: Circularly polarized beam-target double spin asymmetries as a function of 4> at £ = 0.2, —t — 0.4 GeV 2 and Q /2 = 7 
GeV 2 . Left column: Aq x , central column: right column: A© 2 . In the top row panels, 0 integrated over [J, ©p] and the 

calculations are done for different GPD contributions to the TCS process. In the bottom row panels, the calculations are done 
for different 0 angles: 0 — 90° (red dashed curve), 6 — 45° (blue dotted curve) and 0 integrated over [^, (black solid curve) 
and with the GPD H (Reggeized) contribution only. 




Figure 17: Circularly polarized beam-target double spin asymmetries as a function of t. Left: ^4©a;, center: A &y , right: Aq z . 
Calculations are done at the kinematics: £ = 0.2, and Q /2 = 7 GeV 2 and 0=0° (A© x ), 0 = 90° (Aq v ) and <\) — 0° (Aq z ). 6 is 
integrated over [J, in these calculations. 
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Figure 18: Linearly polarized beam-target double spin asymmetries as a function of 0 at £ = 0 . 2 , —t = 0.4 GeV 2 and Q' 2 = 7 
GeV 2 and 6 integrated over [^, ^-] for different GPD contributions. Left column: Ag x , central column: Ag y , right column: 
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Figure 19: Q /2 -dependence (left panel) and |£|-dependance 
(right panel) of the ratio of the 2-fold cross differential cross 
sections da/dtdQ' 2 of Bethe-Heitler and of TCS (integrated 
over the decay angles: (j) E [0, 2tt] and 0 E [ J, ^f]) calcu¬ 
lated with the gauge invariance restoration term (dot-dashed 
curve for TCS: gauge invariant/asymptotic) and with the ex¬ 
act kinematics (dotted curve for TCS: ”exact”/asympototic) 
to the asymptotic limit (i.e. Bjorken limit) calculation. The 
calculation has been carried out for £ = 0.1, 0 is integrated 
over [7r/4, 37r/4], 0 is integrated over [0, 27r] and at —t = 0.4 
GeV 2 (left panel) or at Q' 2 — 7 GeV 2 (right panel). 
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Figure 20: Difference of the single and double spin asymme¬ 
tries (circularly polarized photon) as a function of |t| before 
and after the corrections. Left panel: with the exact kine¬ 
matic compared to the asymptotic limit. Right panel: with 
the gauge invariance restoration term compared to the asymp¬ 
totic limit. The calculations have been carried out for £ = 0.2, 
—t = 0.4 GeV 2 , Q' 2 = 7 GeV 2 and 0 has been integrated over 
[7r/4, 37 t/ 4]. The different (j) angles are indicated on the figure. 















